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Analysis of correlation of multiplicities between various rapidity bins is carried out in the frame¬ 
work of a superposition approach consisting of three phases of the ultra-relativistic nuclear collision: 
early partonic phase, intermediate collective evolution, and statistical hadronization. Simple rela¬ 
tions between the moments of produced hadrons and the moments of the sources produced in the 
initial partonic phase are presented. They involve only a few effective parameters describing the 
microscopic dynamics of the system. We illustrate the practicality of the approach with the Glauber 
model simulations. Our study bears direct relevance for the interpretation of the upcoming results 
for the multibin multiplicity correlations from the LHC, which will help in an assessment of the 
correlation features of the state formed in the earliest stage of the reaction. 
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I. INTRODUCTION 

Long-range rapidity correlations between multiplicities 
of produced hadrons have been under active considera¬ 
tion since the early experiments with pp and pp colli¬ 
sions m, followed with heavy-ion data in the RHIC 
era Numerous theoretical investigations followed, 

in particular [HH22]. The interest is driven by the expec¬ 
tation that such correlations may provide understand¬ 
ing of the elementary space-time dynamics in the earliest 
(partonic) stages of the reaction. An extension of the con¬ 
ventional forward-backward analysis based on multibin 
correlations in rapidity and applying factorial moments 
has been developed in [2SJ [25]. 

We recall that according to the interpretation of the 
STAR Collaboration data of Ref. |2J El , made in [19) 22] , 
we are quite far from understanding theoretically in 
a complete way the forward-backward correlations in 
heavy-ion collisions, in particular their centrality depen¬ 
dence. With the analyses from the Large Hadron Collider 
(LHC) expected to be released in the near future, it is 
important to prepare theoretical ground for the interpre¬ 
tation of the data, which is one of the goals of our work. 

In this paper, which extends the analysis of Ref. [3D] 
to the case of three bins (forward, central, and back¬ 
ward), we analyze the pseudorapidity correlations in a 
superposition approach. In this framework, to describe 
the statistical aspects of relativistic heavy-ion dynamics 
we use a common picture involving three separate stages: 

1) early generation of initial partonic densities (sources), 

2) collective (hydrodynamic or transport) evolution, and 

3) statistical hadronization. A crucial assumption for the 
method is that the emission from a source is independent 
from other sources and universal, i.e., occurs from the 
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same statistical distribution independently of the rapid¬ 
ity or the centrality class of the collision. This leads to 
very simple formulas for the statistical measures, relat¬ 
ing the moments of the initial sources to the moments of 
the distribution of the produced hadrons, in particular 
the multibin correlation coefficients. Parameters present 
in these relations depend on the features of the overlaid 
statistical distributions and properties of hydrodynamics, 
and have simple interpretation. 

We use the derived relations to obtain predictions for 
multiplicities, variances, and multibin correlations of the 
produced hadrons. To model the initial phase, we use the 
Glauber framework based on wounded nucleons |51 [ - 13 31 
amended with binary collisions [31H3H1- The predictions 
are made for the Pb+Pb collisions at y/s^N = 2.76 TeV, 
as the corresponding analyses from the LHC are expected 
to be published in the near future. 


We note that simulations of the forward-backward mul¬ 
tiplicity correlations for the Au+Au collisions at the 
highest RHIC energies were carried out by Konchakovski 
et al. HZ] in the framework of the wounded-nucleon model 
with overlaid Poisson distribution and in the Hadron- 
String-Dynamics transport approach [37]. That work 
carefully discusses the influence of centrality selection on 
the correlation measures used in the experiment, point¬ 
ing out to effects stemming just from the fact that the 
centrality classes are wide. Our results from the Glauber 
model comply to these general findings. 


Our paper is organized as follows. In Sec. [TT] we re¬ 
view the three-stage approach, consisting of the early 
production, the hydrodynamic evolution, and the final 
freeze-out. We derive the relations between the multi¬ 
bin moments of the sources produced in the early phase 
and of the multibin moments of final hadrons. Some 


general predictions of the approach are given in Sec. Ill 


In Sec. [TV] we discuss the technically important issues 
related to the choice of the centrality bins. For the con¬ 
sidered statistical measures, scaling laws with the width 
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FIG. 1. (color online) Time evolution of the created system, 
visualizing the subsequent stages of particle production. B, 
C, and F indicate the backward, central, and forward bins, 
respectively, extending along the spatial rapidity or pseudo¬ 
rapidity direction. 



FIG. 2. (color online) Initial source projected on the trans¬ 
verse plane (left) and the three dimensional picture of color 
flux-tubes extending along the spatial rapidity (right). 


In our considerations it is only important that h is 
(to a good approximation) proportional to p. 

3. Freeze-out and the statistical hadronization mech¬ 
anism (see, e.g., |47| for a review), converting 
the hydrodynamic fluid into hadrons n, which are 
then subject to experimental observation in various 
pseudorapidity bins. 


of the centrality bin are derived. Section |V] presents our 
results of the Glauber simulations for the initial phase. 
We show here the practicality of our method in the appli¬ 
cation to the future experimental data, which will serve 
to verify the model of the initial phase. 


II. THREE STAGE APPROACH 


The above simple scheme summarizes most of the pop¬ 
ular model approaches to particle production in rela¬ 
tivistic heavy-ion collisions. Our treatment focuses on 
statistical aspects of each of the phases. While the un¬ 
derlining dynamical mechanisms are very complicated, it 
turns out that certain statistical features may be simply 
parametrized, and specific conclusions concerning parti¬ 
cle correlations in rapidity may be, under mild assump¬ 
tions, drawn in a model-independent way. 


In Ref. [HQ] we have considered a three-stage approach, 
based on superposition, that describes certain features of 
production of particles in relativistic heavy ion collisions 
and their correlations, focusing entirely on statistical as¬ 
pects. In this framework, the complicated dynamical 
evolution of the system is represented with the follow¬ 
ing separate stages: 


1. Initial phase, which may be modeled with the 
Glauber approach [3lH34l l36l 138] , color glass con¬ 
densate [5D1II2] , string formation mm, etc. As the 
result of this early dynamics, some initial space- 
time distribution of entropy or energy-density is 
generated. We assume that first the initial sources 
are formed (see Sec. |II A ), denoted with s, which 
then produce the distribution of initial partons, p 
(see top of Fig. [l]). The extent of correlation be¬ 
tween the source densities in spatial rapidity bins 
B , C , or F depends on the early dynamics. 


2. Hydrodynamic evolution (for reviews see, e.g., [441 
I45| and references therein) or transport [461 evolves 
the initial density of partons p into expanding hy¬ 
drodynamic fluid density h (third row in Fig. [I]). 


A. Sources 

We use the concept of the initial sources , which yield 
the initial space-time entropy or energy density distri¬ 
bution. The precise definition depends on the model 
of the initial phase. In our work we use the Glauber 
model [3IH33] ESI HU US], where the source is created 
in an individual NN collision, located in the transverse 
plane and extending along the spatial rapidity direc¬ 
tion, 77|| = \ log[(t + z)/(t — z)\ (cf. Fig. [ 2 ]). A particular 
model of the distribution in r/y, coming from Ref. [3D] and 
used here, is described in Appendix [A] It is motivated 
with the success of Ref. m in describing the deuteron- 
gold collisions [35], followed with the extension to NN 
collisions US |3S [51] . We note that other descriptions of 
the early production could be used as well, for instance 
the Kharzeev-Levin-Nardi framework [34l 1331 [56] , 

In certain models, the sources may be interpreted as 
color flux-tubes, extending along some range in rapidity, 
as pictured in Fig. [ 2 ]). It may happen that some flux- 
tubes are longer, extending through the whole region of 
interest from the B to F bins, while some may be shorter, 
as indicated in Fig. [2j) , making the correlations shorter- 
range. Clearly, such effects influence the multiplicity cor- 














































3 


relations between the bins, which is what ultimately we 
wish to investigate. 

The dots in Fig. [2] indicate the sources in bin A, where 
A = F, B, C. Their corresponding numbers are denoted 
as sa throughout this paper. 

B. Centrality definition 

The number of sources in the central bin, sc, is used 
to define the centrality classes of the collision in our sim¬ 
ulations presented in Sec. [V] In the wounded nucleon 
model this prescription would correspond to using the 
number of wounded nucleons to determine centrality, as 
is frequently done in model studies. The definition based 
on wounded nucleons or participants of the projectile is 
appropriate for experiments which measure the projec¬ 
tile spectators in dedicated detectors (e.g, as in NA61). 
On the other hand, in experiments where the multiplic¬ 
ity of the event is used as the centrality selector, a closer 
measure would be the number of hadrons generated in 
the model, e.g., nc■ As pointed out in Ref. m the 
definition of centrality may influence the results for the 
forward-backward correlations to some extent. We inves¬ 
tigate this issue in Sec. |V C| 

C. Initial dynamics 

In the early stage of the reaction, the sources emit 
partonsQ We assume that each source produces partons 
independently of one another and, in addition, the pro¬ 
duction occurs independently in the F, B , and C bins, as 
they are assumed to be well separated in the spatial ra¬ 
pidity 77 || . Thus, each source i emits pi partons according 
to the same distribution, 

SA 

Pa = £>, A = B,F,C. (1) 

2=1 

It is then elementary to obtain the formulas for the above 
superposition model for moments corresponding to aver¬ 
ages over events (here listed up to rank three): 

(pa) = ( p)(SA), 

( A Pa) = var (p)(s A ) + (p) 1 2 (As 2 A ), 
(Ap A Ap A ') = (p) 2 (As a As a'), 

( Ap \) = Paip) (sa) + 3var(/i)(/z)(A s\) 

+ (m) 3 (a4), (2) 

(Ap A Ap A ,) = var(p)(p)(As A AsA>) 

+ ( p) 3 (As a As 2 a ,), 

(Ap F Ap c Ap B ) = (p) 3 (As F As c As B ), 


1 We note that this mechanism was implemented already in the 

original wounded nucleon model m , where the Poisson distribu¬ 

tion was overlaid over the distribution of the wounded nucleons. 


where A, A' = F, £?, or C, and A ^ A! . In our notation 
A p A = (p A — (p/ i)) n , n = 1,2,3, while (•) denotes aver¬ 
aging over events. The symbol p?,(x) stands for the third 
central moment of the distribution of x. 


D. Hydrodynamics 

The density of partons formed in the early phase, as 
described above, provides initial conditions for the equa¬ 
tions of the hydrodynamic evolution of the fluid cell in 
the intermediate phase. This evolution, proceeding until 
freeze-out, is assumed to be deterministic (note, however, 
Refs. ISZllSH]), hence no extra fluctuations are introduced 
at this stage. As a result, the initial density of partons p 
is deterministically carried over to the entropy density of 
the fluid cell at freeze-out, denoted with h (cf. Kig. [l]). 

In the vicinity of the average value of p the function 
h(p) is approximately affine, hence we may expand 

h = t 0 (p)+ti(p- (p)) + 0((p- (p)) 2 ) , (3) 

where U are parameters depending on the dynamical re¬ 
sponse of hydrodynamics to the initial condition. The 
higher-order terms in Eq. § may be dropped if p is suf¬ 
ficiently close to (p), otherwise corrections arise to the 
formulas developed below. With Eq. ([3]) we get 

(h A ) = t 0 (pA), 

(.Ah 2 A )=t\(Ap \), 

(Ah A Ah A >) = t\(Ap A Ap A >), 

(Ah\) =t\(Ap\), (4) 

(Ah A Ah A ,) = tl(Ap A Ap A ,), 
(Ah F AhcAh B ) = t\(Ap F ApcAp B ). 

We note that the approximate linearity of the hydro- 
dynamic response has been well determined for the evo¬ 
lution of eccentricities [591 ED] at a given centrality. For 
the present case of multiplicities, the success of modeling 
of the measured particle multiplicity distributions within 
the Glauber framework [Ml ED indirectly supports the 
assumptions leading to Eq. Q. 


E. Statistical hadronization 

When hydrodynamic evolution ends, statistical 
hadronization phase takes over at freeze-out (cf. Fig. [TJ. 
We assume that the well-separated fluid cells h emit n 
hadrons independently from one other into correspond¬ 
ing well-separated phase-space regions in pseudorapidity. 
Each of the h cells emits m hadrons in the same universal 
manner, hence 

hA 

n A =^2, m i- (5) 

2=1 
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Then, in full analogy to Eq. ([2]), we obtain 
(n A ) = ( m)(h A ), 

(A n A ) 2 = var (m){h A ) + (m) 2 (Ah A ) 2 , 
(An A An A >) = ( m) 2 (Ah A Ah A > }, (6) 

(An\) = ^3 (m) + 3var(m)(m)(A/r^) 

+ (m) 3 (Ah A ), 

(An A An%) = var(m)(m) (Ah A Ah A >) 

+ (m ) 3 (Ah A Ah A ,), 

(AnFAncAns) = (m) 3 (A/i^ AhcAhs)- 

The numbers nn are observed experimentally. We note 
that the distribution of mi incorporates the effects of the 
detector acceptance. 


The scaled variance, used in the following, is defined as 
u(x A ) =var(x A )/(x A ). (12) 

We also introduce the combinations of parameters 

1 2 

6 = /3/a = oj(m) H— ~(m)u>(p), 
to 

t 2 

K = 1 /a= ±(n)(m), (13) 

to 

t 0 u](m) w(p) 

III. GENERAL RESULTS 


F. Combined formulas 


We are now ready to combine Eqs. PPP , which yield 
relations between moments of produced hadron distribu¬ 
tions and moments of the initial sources: 


(n A ) = a(s A ), 

(A n\) = /3(s A )+"f(As 2 A ) 

(An A An A ') = 7(AsaAsa'), (7) 

(A n\) = C(s A ) + 3j 1 / 2 p'(A s\) + 7 3/2 (As^), 
(An A An 2 A ,) = ^ 1/2 /3’(As a As A ') +'y 3 / 2 (As A As 2 A ,), 
(AnpAncAiiB) = ^^(AsfAscAsb)- 

The constants introduced above are defined as 


( 8 ) 


a = t 0 (p)(m), 

/ 3 = fo(/t)var(m) + t 2 (m) 2 v&r(p), 

(3 1 = ti(/x)var(m) + f 2 (m) 2 var(/x), 

1 = tl(m) 2 (ii) 2 , 

C = toii 3 (m)(n) + 3t 2 var(^)var(m)(m) + tlp 3 (p)(m) 3 . 

The algebraic structure of Eqs. 0 is very simple. How¬ 
ever, nontrivial character stems from the fact that the 
parameters <[8j) do not depend on centrality of the colli¬ 
sion. 

It is convenient to introduce a short-hand notation 


Sijk — (As x b AsqAsf) 
for the moments of rank r = * + j + k : and 


Sijk — 


(As l B As J c Asp) 

V(A4)(A4)(A s ly 


(9) 


( 10 ) 


for the moments scaled with the standard deviations. In 
particular, the forward-backward correlation of the num¬ 
ber of sources is 

ri _ (AsbAsf/ _ . . \ 

Jioi = ~ 7 = . , = p(SF, SB)- (11) 


V( As 2 b)(As 


A. Relations of moments of hadrons to moments of 
sources 

Our methodology uses the Glauber model for the 
sources to make predictions for the moments of n A . 

An interesting relation for the hadron correlations re¬ 
sults from Eq. ([7]): 


p(n A ,n A >) = p(s A ,s, 4')l 


1 - 


4»a) 


1 - 


j(n A ') 


(14) 


Since, as we shall see, for not too large pseudorapidity 
separations, in the Glauber model to a high accuracy 
p(s A , s A ') ~ 1, it then follows that 


p(n A ,n A >) ~ 


1 - 


u(n A ) 


1 - 


u)(n A ’ 


■ (15) 


which means a very specific relation between the correla¬ 
tion of numbers of hadrons in bins A and A' and the cor¬ 
responding scaled variances, which depends on a single 
centrality independent parameter <5. Relation (15) may 
be straightforwardly tested with the future data. 

From Eqs. ([TJ) we also find 


uj(n A ) = 5 + ku>(s a ), 

p(sf, Sb ) 


p(n F ,n B ) = 


1 + 


A 

o(sf) 


(16) 


1 + 




formulas presented already in Ref. [[30 1. 

Another interesting relation, linking three bins, has the 
form 


(AnsAncAnF) (AsbAscAsf ) 
(An A An A >) 3 / 2 (As a As A ') 3 / 2 


(17) 


Again, the right-hand side, obtainable in a model, can 
be verified against the future data, yielding the left-hand 
side. 
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With the help of the basic formulas (J7| it is straight¬ 
forward to find relations for the newly-proposed mea¬ 
sures |B2, .53] of fluctuations. They are defined as 


S (n F ,n B ) 


A (n F ,n B ) 


('n F )uj(n B ) + (n B )uj(n F ) - 2 co v(n F ,n B ) 


Cs 

(n B )u(n F ) - ( n F )uj(n B ) 

C A 


(18) 


where the constants Cs and Ca are proportional to first 
moments of any extensive quantity. We take the simplest 
choice 


Cs = Ca = (n F ) + (n B ). (19) 

Elementary algebra leads to the relations 


E(riF, n B ) = S + kE(sf, s b ), 

A(n F ,n B ) = ^ +kA(s f ,s b ), ( 20 ) 

( s F + s B ) 


where T,(n F , n B ) and A (n F ,n B ) are defined analogously 
to Eq. (18), but with ua replaced with s^. 


B. Implications for the qualitative features of 
hydrodynamics 


Based on first and third of Eqs. ([7|, we find the relation 

to 
ti 


7 


a 

~T/2- 


( 21 ) 


Since the a and 7 parameters may be found by fitting 
the model predictions to the future data, we are able to 
obtain from Eq. (21) information concerning the hydro- 
dynamic response in the intermediate phase. Figure ([3]) 
shows two possible scenarios for the hydro evolution. 
When to > t\ (to < ti), the hydrodynamic growth h(p) 
is slower (faster) than the linear function. From the ex¬ 
perience of hydrodynamics we expect that to — ti, which 
means that the entropy in the considered cell is propor¬ 
tional to its initial entropy (see, e.g., IMll. 


IV. DEPENDENCE OF RESULTS ON WIDTH 
OF CENTRALITY BINS 

An important question concerns the dependence of the 
results shown in this paper on centrality selection. In 
our analysis presented below we define centrality with 
the help of the number of sources in the central bin, sc- 
The results for moments of multiplicities depend on the 
chosen width of the centrality bin Ac - in that sense 
they are not intensive measures [62( I65|. Nevertheless, 
the moments evaluated that way are still useful for com¬ 
paring model predictions to experiment, as exactly the 
same choice of width of the centrality windows can be 
made in model simulations as in the data analysis. 

To understand the problem, let us focus on Fig. [4j 
where we show a result of a simulation in the form 



FIG. 3. (color online) Possible scenarios for the hydrodynamic 
response (see text for details). 


of a scattered plot of s F and s B . The short skewed 
lines cut the events into subsequent centrality classes 
with Ac = 2.5%, while the longer lines into classes with 
Ac = 10% (containing 4 bins of Ac = 2.5%). We note 
that the scattered plot is very much elongated, sim¬ 
ply reflecting the large correlation between s F and s B . 
With the exception of very small centralities c (the lower 
left corner of the plot), the correlation between s F and 
s B within a given centrality window is very close to 1 , 
p(s F ,s B ) ~ 1. As we move towards very low c (lower- 
left corner), or when Ac is decreased, p(s F , s B ) becomes 
lower than 1. Thus the correlation of the sources is in an 
obvious way affected by the choice of the centrality bins. 

For sufficiently large width of the bins, the above- 
discussed feature allows us for a derivation of approxi¬ 
mate scaling laws of the moments on Ac. Details are 
shown in Appendix [A] The average multiplicities do not 
scale, as long as the window is sufficiently narrow: 

(a a) ~ (Ac) 0 . ( 22 ) 

For the variance and covariance we find 

((As a ) 2 ) ~ (Ac) 2 , 

(As A As A ') ~ (Ac) 2 , (23) 

from where 

p(s A ,s A ') ~ (Ac) 0 . (24) 

For the moments of rank 3 we need to be a bit more 
sophisticated. The marginal distribution for the variable 
sa in a given centrality bin can be approximated as an 
affine function, /(sa) = N(1 + as a), where a is small. 
Then (cf. Appendix [A|, to leading order in a 

((Asa) 3 ) ~ «(Ac) 4 , (25) 

and 

((As^) 3 )/ (As a Asa') 3 / 2 ~ aAc. (26) 

The above scaling laws are visible in our results presented 
below. 
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FIG. 4. (color online) Scattered plot for the forward-backward 
correlation of the numbers of initial sources sf and sb for 
rapidity separation A rj = 1.8. The short lines cut the sam¬ 
ple into centrality bins of width Ac = 2.5%, and the long 
lines into bins of width Ac = 10%. GLISSANDO, mixed model, 
Pb+Pb collisions at the LHC energy of y/ snn = 2.76 TeV. 


V. GLISSANDO SIMULATIONS 

We use GLISSANDO [BB] to model the initial 
phase. To simulate Pb+Pb collisions at the LHC at 
yjS nn = 2.76 TeV, we use (unless stated otherwise) the 
mixed Glauber model j34H36] with the mixing param¬ 
eter a = 0.1 (cf. Appendix [5]), the NN cross section 
&nn = 66.3 mb and a Gaussian wounding profile EH- 
The numbers of sources sa are obtained in pseudo¬ 
rapidity bins of width Srj = 0.2, while Ap denotes the 
separation of the centers of the F and B bins. The C 
bin is centered around zero, while the F and B bins are 
arranged symmetrically. 


A. Results for the correlations of sources 


We begin presenting our results with the forward- 
backward and forward-central correlation functions, 5ioi 
and Sno, shown in Fig. [ 5 ] for the rapidity separations 
A 77 = 1.8 and A// = 10. We use Ac = 10 %, but as 
follows from Eq. (23), the results are insensitive to the 


width of the bin. We note that these correlation func¬ 
tions are essentially equal to unity, except for very large 
values of A 77 and c. This is an interesting feature of 
the adopted Glauber model, as it shows that we can use 
p(sa,sa') — 1 in practical applications. This fact has 
its origin in a strong correlation between the numbers 
of wounded nucleons from nuclei A and B, as well as be¬ 
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FIG. 5. (color online) Forward-backward and forward-central 
correlations of the number of initial sources, plotted as func¬ 
tions of centrality for two cases of the pseudorapidity separa¬ 
tion, Arj = 1.8 and A 77 = 10. 



FIG. 6 . (color online) Scaled rank-3 moments of sources plot¬ 
ted as functions of centrality for three sample values of the 
width for the centrality bins. 
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FIG. 7. (color online) Prediction for the combination of 
third and second moments of observed hadrons, plotted as 
a function of centrality for three pseudo-rapidity separations, 
At] = 1.8, 6 , 10. 


FIG. 8 . (color online) Same as in Fig. [t] but for different 
variants of the Glauber model. 


tween the wounded nucleons and the binary collisions (cf. 
Appendix [B]). 

For the case of symmetric collisions and symmetri¬ 
cally placed F and B bins, the measures Y,(sf,Sb) and 
A (sf,sb) from Eq. (18) are equal to 


£(sf, s b) = w(s^)[l - p(sf, sb)], (27) 

A(sf, s_b) = 0, (A = F,B), 

Since p(sf,Sb ) — 1, it follows that T,(sf,sb) — 0. 

The scaled rank-3 correlations of Eq. (10) are displayed 
in Fig. [ 6 ] We note the advocated scaling with the value 
of Ac, cf. Eq. (26). Moreover, various rank-3 moments 


shown in the plot are essentially equal to one another in 
the Glauber treatment. 



B. Results for the correlations of produced hadrons 

Now we pass to the presentation of results for the mo¬ 
ments of the distributions of the produced hadrons, tia- 
Under the assumptions of the superposition approach, 
Eq. ( p~7] ) allows us to make predictions for the produced 
hadrons, based on the model of sources. Our result, 
for several pseudo-rapidity separations Ar/, is shown in 
Fig .0 The result is specific to the adopted source pro¬ 
duction model. In Fig. [ 8 ] we compare the outcome, for 
the same quantity, of the wounded nucleon model (mix¬ 
ing parameter a = 0 ), the mixed model (a = 0 . 1 ), and 
the binary-collisions model (a = 1). The large differences 
between the models should allow, with the future data, 
to discriminate between them. 

Finally, we use Eq. (J7| with some reasonable choice of 
the parameters d8|T3] ) to make predictions for the multi¬ 
plicity, variance, covariance, and rank -2 correlation coef¬ 
ficients. The results are displayed in Figs. m We may 


FIG. 9. (color online) Predictions for the multiplicity of the 
produced hadrons in peripheral bins, plotted as a functions of 
centrality, based on Eq. ffl at different values of parameters 
in Eq. ppt - 


now carefully illustrate our methodology: with the data 
for multiplicity we fit the value of a with the results of 
Fig-! the value of 7 with Fig. [l2j and then the value of 


f$ with Fig. [IOJ Next, we evaluate A = /3/q and make a 
test with the help of Fig. [13] Extension to rank-3 mo¬ 
ments is straightforward. That way, we may obtain the 
parameters of the superposition model in a systematic 
way. Knowledge of these parameters will provide insight 
into the microscopic features of the particle production 
mechanism in the system. 
























FIG. 10. (color online) Same as Fig. [9j but for the variance. 



FIG. 11. (color online) Same as Fig. [ 9 ] but for the scaled 
variance. 



FIG. 12. (color online) Same as Fig. [ 9 ] but for the covariance. 



FIG. 13. (color online) Same as Fig. [9] but for the forward- 
backward correlations. 


C. Centrality from the number of hadrons 


In this Section we check how the results for the 
forward-backward correlations depend on the method of 
determining centrality, carrying out a study similar to 
Ref. |T7j. In Sec. IIB we have introduced the definition 
of centrality, used in the remaining parts of this paper, 
based on the number of sources in the central rapidity 
bin, sc- To estimate the results based on centrality de¬ 
termined with number of produced hadrons in the central 
bin, nc, we carry out GLISSANDO simulations, where a 
negative binomial distribution is superimposed over the 
distribution of sources. The parameters of the binomial 
distributions are such that the experimental multiplicity 
distributions at the LHC are reproduced [ 68 ] , 

The qualitative effect of the superimposed fluctuations 
is simple to understand: when the limits of the central¬ 


ity bin are n}™ and n^ lgh , the corresponding limits in the 
number of sources, s^ w and Sp lsh are somewhat diffused. 
This causes more extended scattered plot of type of Fig. [4] 
for the events belonging to the considered centrality class 
and the correlation coefficient p(sf,Sb) increases. As 
may be inferred form Fig. |4j the effect will be larger for 
peripheral collisions, where the scattered plots are less 
elongated, and smaller for central collisions, where elon¬ 
gation (for Ac = 10%) is large. 

Results of our numerical simulations are presented in 
Fig. [l4j where we show the ratio of the correlation co¬ 
efficients 5ioi for the centralities evaluated according to 
sc and according to nc obtained with the overlaid neg¬ 
ative binomial distribution, indicated as Sioi/S^. We 
note that for small forward-backward pseudorapidity sep¬ 
aration, A= 1 . 8 , there is no visible effect, whence 
S 101 /S 101 — 1. For large A 77 and peripheral collisions 
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FIG. 14. (color online) Influence of the method of central¬ 
ity selection for the forward-backward correlation of sources, 
tested with the ratio 1 S 101 /Sioi of the forward-backward cor¬ 
relations evaluated with the number of sources sc and the 
number of hadron nc in the central rapidity bin. At small 
forward-backward separation, A?; = 1 . 8 , there is no effect of 
the method used to determine centrality. 


S-f oi /•S'un < 1 , in agreement with the qualitative argu¬ 
ments presented above. 


VI. CONCLUSIONS 
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Appendix A: Centrality width dependence 


In this Appendix we sketch the derivation of scaling of 
the moments with the width of the centrality bin. The 
explanation comes from Fig. 0 , where the forward- 
backward correlation is plotted as scattered plot. For 
a given centrality and for a typical size of Ac we may 
approximate the distribution of = x as an affine func¬ 
tion, 


fix) 


1 + ax 
(1 + ac ) Ac’ 


(Al) 


which has been normalize to 1. Upon expanding for small 
ac it follows that 

r c+ ^r a Ac 2 

(x)= x/(x)^c+—, (A2) 


((Ax) 2 ) = 


rc+^f 


(x-(z)) 2 ^^, 

Ac 1 Z 


((Ax) 3 ) = 




(x - (x)f ~ Ac 4 , 


((Axf) 

((Ax ) 2 ) 3 / 2 


Ac. 


The scaling laws (|22|23 ) result from the above formulas. 


We have presented a systematic study of two- and 
three-rapidity-bin multiplicity correlations in relativistic 
heavy-ion collisions in a superposition approach, based 
on three phases of the collision: initial partonic phase, 
intermediate collective evolution, and final statistical 
hadronization. The derived expressions show how the 
moments of the produced particles are related to the 
moments of the sources produced in the initial phase. 
The relations, simple but nontrivial, involve only a few 
parameters collecting the microscopic information con¬ 
cerning the dynamics (parameters of the superimposed 
distributions, response of hydrodynamics). 

With a definite model for the production of sources, 
such as the Glauber model used here, one is then able to 
obtain in a systematic way predictions for the multibin 
multiplicity correlations of the produced hadrons. Con¬ 
versely, with the future experimental data, one will be 
able to test the models of the initial phase, as well as ob¬ 
tain information on the microscopic parameters. In par¬ 
ticular, the awaited forward-backward multiplicity cor¬ 
relation analysis with the LHC data will shed light on 
the early production mechanism. The statistical method 
presented in this paper is directly applicable to that case. 


Appendix B: Emission profiles in spatial rapidity 


In a nucleus-nucleus collision, the emission profile [50] 
defining the shape of initial sources in the spatial rapidity 
has the form 



(/+ fa |\) N w , b + f-(v\\) n -w,f) 


+ a 


^ /+fall) + /-fa|| 


N, 


bint 


(Bl) 


where N w p and N w p are the numbers of the forward- 
and backward-going wounded nucleons, Nf )in is the num¬ 
ber of binary collisions, while /+ (? 7 ||) and /_ (? 7 ||) de¬ 
scribe the corresponding wounded-nucleon emission pro¬ 
files. Finally, / ( 7711 ) is the emission profile for the binary 
collisions. We use the following parametrization EHj: 

/ fall) = ex P^— 6>(|77||| - r?o)^ , 

/+fall) = .frfall) / fall) > (B2) 

/- fan) = fF (~V\\) f (v\\), 
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FIG. 15. (color online) The emission profiles in space-time rapidity for the wounded nucleons (dashed lines) and the binary 
collisions (solid line). The profile /+(/_) corresponds to the forward (backward) moving wounded nucleons. 


with 


fp fan) 


0 , rm - 

-Vm < V\\ < Vm 

1 r)m<V II- 


T]m = 8.58, a n = 1.4. 
Fig. [15} 


The functions (B21 are shown in 


Sample values of parameters [55], describing the ATLAS 
data after the hydrodynamic evolution, are v o = 1-25, 
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